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1 Introduction 



In the last few years two-dimensional Yang-Mills theories have been studied exten- 
sively. The partition function for the theory on E fl , a Riemann surface of genus g, was 
computed exactly in |1] (see also 0). More recently Gross and Taylor have provided 
a string interpretation of the 1/N expansion of the partition function for the SU(N) 
gauge group by enumerating classes of maps from worldsheets to the target surface 
|> HI- The string description of the theory was later extended to the SO(N) 
and Sp(2N) gauge groups in ||. Some progress has been made towards finding a 
string action for the 2d QCD string |7], || as well. However Douglas and Kazakov 
|T9| have shown that the string interpretation no longer holds for weak coupling when 
S 9 is a sphere, although it does remain valid for strong coupling. On the sphere the 
large N theory undergoes a phase transition at a critical value of e 2 (Area) = it 2 . 
The study of this phase transition might also be relevant for four dimensional QCD, 
which may or not exhibit a similar behaviour []19, ZOfl . It is therefore important to 



understand those features of a gauge theory which give rise to string behavior and to 
a phase-transition. 

Two dimensional Yang-Mills theories have also been studied by means of beau- 
tiful path integral methods. In || the path integral was evaluated by a generalized 
localization formula, while in |T(| , the path integral for gauge fields of fixed holonomy 



around the boundary of a disc was used as the elementary building block from which 
the partition function and Wilson loops on closed surfaces could be constructed. 

In this paper we will study correlation functions of field strengths in 2d non- 
abelian gauge theories as an application of the abelianization technique for path inte- 
grals developed by Blau and Thompson |ll], [12], [13| . We show that these correlators 
exhibit the almost topological nature of the theory, and can be expressed in terms of 
the higher order Casimir operators of the gauge group. They can also be related to 
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the generalized 2d QCD of [18, 25 1 where higher order Casimir operators appear in the 
partition function. In section 2 we use abelianization to compute the partition func- 
tion in the presence of an external source. This allows us to calculate the correlators 
of field strengths, which is presented in section 3. One might expect the correlators 
of field strengths to be trivial, since the Lagrangian with a source coupled to F^ v is 
quadratic in the fields. However, because of the topological obstructions encountered 
in the abelianization, the correlators are highly non-trivial. In section 4 we show that, 
on the sphere, the regularized correlators of arbitrary numbers of field strength oper- 
ators have second-order phase-transitions in the large N limit. Section 5 shows how 
contractible Wilson loops can be computed from our correlators, in agreement with 
the known formulae for Wilson loops. Our results suggest that, in the abelianization 
gauge, the usual Stokes' theorem can be used to relate the Wilson loop to the field 
strength correlation functions when the gauge group is SU{2) [see equation fl53|)]. In 
Appendix A we compute correlators on the disc in a different gauge, while Appendix 
B is devoted to some Lie algebra details. 

The study of the correlators of the field strengths in the abelianization gauge, 
enables us to calculate the master field for the regularized field strength in the 
large N limit on the sphere, in a sequel to this paper fl3"2"|. Since, as we show in 



this paper, in this gauge these correlation functions are essentially topological, i.e. 
independent of position on the Riemann surface, the master field will also have this 
property. This enables us to construct the master field for itself on the sphere 



| 32| . In fact, the master field which exhibits the whole structure of the (unregulated) 
correlation functions, as is needed for the Wilson loop, can be obtained in a Hilbert 
space representation p6|. 



One hopes to be able to calculate fermion correlation functions on the sphere 
in the large N limit, using the master field for A^. Thus, this paper serves as the 
first step in this program of understanding the coupling of fermions to 2D Yang-Mills 
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theory on Riemann surfaces. To date this has only been acomplished on the plane 
P5| . We are presently studying this application of our results. 

2 The Path Integral Z^ g (J) 

In this section we consider the field strength of the U(N) gauge theory on the compact 
surface of genus g,T, g , coupled to an external source J(x) transforming in the adjoint 
of the gauge group. The path integral describing this situation will be evaluated using 
the elegant abelianization method of ET], O. The resulting partition function Zj] g (J) 
will then enable us to compute the electric field correlators in the next section. 

The partition function in question is, 

Z* g {J) = J VA,eM^ 2 k g ^Tr(£ 2 ) + J Sg d^Tr(J^)] (1) 

where the scalar fields £ a (x) are defined by F^i^x) = (,(x)^J g(:r)e^, where £ = T a £ a , 
with T a a generator of U(N), d\i = ^Jg(x)d 2 x a Riemannian measure on £ 9 , is 
the usual antisymmetric tensor with e i = 1, and \J g(x) is the square root of the 
determinant of the metric on E g . 

The starting point is the first order form of the path integral, 

2 

Zx g (J) = I V^VA^ V \-\k 3 ^Tr(0 - iJf - if Eg d/iTr(0£)] (2) 

where is a Lie algebra(^) valued scalar field also transforming in the adjoint under 
the gauge group. Performing the gaussian integration over <p gives back (|l]). We 
decompose Q into, 

G = t®k (3) 

where t = h@ u(l) generates the usual maximal torus (T) of U (N), with h being the 
su(N) Cartan subalgebra (su(N) denotes the Lie algebra of SU(N), etc). 
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The path integral @ then becomes the product of the integrals over the su(N) 
and u(l) valued fields, where the two sectors will be related by a topological selection 
rule as explained aheacf]. We will evaluate each of these separately beginning with the 
£7(1) sector. Bold type style quantities will be su(N) valued and primed quantities 
will be u(l) valued. In order to account for the nontrivial U(N) bundles over £ 9 , we 
use the topological condition for the u(l) component of the gauge curvature 

Je s dug = 27r 77^f' m an integer. (4) 



It is convenient to decompose masm = m-N + p where m is an integer and p G Z^. 
The topological selection rule connecting the u(l) and su(N) sectors will depend on 
p and therefore we keep p fixed for now. The sum over all possible values of p will 
be performed when we examine the su(N) sector. To impose condition (PS) we insert 
the periodic delta function, 

Z) ex P(*-^Js 9 ^0 ■ exp(-27ri— n) (5) 



in the path integral. Let us consider the w(l) piece first, using the strategy developed 
in fLOf . The action is linear in £' , and the Nicolai map can be used to change variables 
of integration from A'^ to £ , and some scalar gauge fixing function G(A' IM ). The 
integration over this second new variable, together with the jacobian from the change 
of variables, will cancel against the Faddeev-Popov determinant. The integration 
over £' now gives a delta function that fixes the (constant) value of 0'. The final 
contribution to (§) from these u(l) abelian fields is, 

£ ex P(^-/s 9 d V ( J 'f - + i -7^ e2n k g J ') ■ exp(-27rz— n) (6) 



where the source J (x) just couples to u(l) fields^, and A is the area of S g . (We 
take Tr (T a T b ) = 5 ab for the U(N) generators). Notice that for N = 1 there is no 

1 This condition reflects the fact that we have U(N) gauge group, not just U(l) x SU(N). 
3We use / = ^ and { = . 
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dependence on p in (^). We will now do the remaining integrations over the su(N) 
valued fields. Here we will use the abelianization method of [O, |T!1. One chooses the 



gauge condition, which restricts <fr to the Lie algebra of the maximal torus, i.e. 

<f) k = (7) 

in order to make use of the fact that the inner product (given by Tr) on Q makes h 
and k orthogonal. First one must cover E 5 by open sets where (0) is valid, then in 
the intersections of those open sets with each other, the h components of <fi (and also 
other fields) will in general be related by nontrivial gauge transformations preserving 
([/]) i.e. by transformations which take values in N(T), the normalizer of the maximal 
torus T in G. The remarkable fact that one can choose those gauge transformations to 



lie in T itself [I3|, so that is globally defined, follows from the simple connectivity 



of Q. These T-valued gauge transformations then give rise to non-trivial T bundles 
over Tig. We refer to Blau and Thompson's [pl| for the detailed explanations of this 



beautiful result. (Strictly speaking at this point one must restrict the fields to take 
values in the regular elements of Q, but we'll see ahead that the non-regular valued 
fields are eliminated by the Faddeev-Popov determinant [13]). Thus, the price to be 
paid for demanding (|7j) is the appearence of nontrivial T bundle topologies. These 
however can be easily incorporated in the path integral by making use of the first 
Chern numbers [[11], |12|]. We have then, for the integration over the su(N) valued 
fields 

L £ JvA^VA^detlAwifyeMSiJ)) (8) 



Tbundles 



h 

where \W\ is the order of the Weyl group and det[A^(0 )] is the path integral 
analogue of the Weyl determinant in the Weyl integral formula for Lie groups. This 
can also be interpreted as a Faddeev-Popov determinant coming from ([7|), ||12|| . 

det[A w (/)] = det[ad(/)| n0(S k) ] (9) 



where f2°(£ 9 , k) is the space of k valued 0-forms on E 
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Using the orthogonality between h and k the only contribution from the inte- 
gration over the k components of the gauge connection is easily seen to be given, up 
to a constant, by 



(detM^)^^])- 1 / 2 (10) 

where f2 1 (S 9 , k) is the space of k valued 1-forms on E g . Notice the similarity to 
In fact using the Hodge decomposition of k) into orthogonal components, the 

contribution of (|9]) and (|i0|) nearly cancels (see for a careful evaluation of this 
ratio of determinants); that is 
h 

<l0t ' [ '" ,{ ° )Ll " '' J detM/)y-^» = detM/)y^)/ 2 (11) 



(det[ad(^)| nl(I!jJb) ])^ 

where 60 — 1, &i = 2g are the Betti numbers for E ff and x(^g) is the Euler number of 
the surface. (We'll see that only constant (jfr configurations contribute in the end). 



At this stage we have reduced the integration over the sw(iV)-valued fields to 



Tbundles ' 



(12) 



The remaining T valued gauge invariance, not eliminated by (|7|), can be fixed, 
and the corresponding Faddeev-Popov determinant eliminated by using the Nicolai 
map once more. In this way dA^ 1 is traded for A^ as a new variable in the path 
integral. It is now easy to express the summation over nontrivial T bundles in terms of 
the first Chern numbers, as in (|J) and (||). We choose a set aij, i — 1, N—l, of simple 
roots for su(N) and we let (dA^ 1 ) 1 be the components of dA^ 1 in the corresponding 
basis for h. Let us consider the topological sector with u(l) charge m as in (|4]). If we 
write m = fh ■ N + p with p e and fh an integer, then the topological conditions 
for (dA^y are 

k g (dA h y = 2ttK - ^), w? integers (13) 
7 



This topological selection rule arises because the U(l) and SU(N) components of 
U(N) are related since U(N) ~ (U(l) x SU(N))/Z N . Inserted in the path integral, 
this condition will combine with the exp(— 2m^n) factor in @ and will give origin 
to a summation over all pairs of representations of U(l) and SU(N) that can be 
combined into a representation of U(N), namely such that (n — R) modiV = 0, where 
R is the number of boxes of the Young tableau for the SU(N) representation and n 
labels the U(l) representation as in (||). Therefore the summation over T bundles 
and the remaining summation over p can be done by inserting 



£ exp[*A(/ E9 dA n )} ■ 6((n - R{\)) modTV) (14) 
AeA 



in the path integral. In the sum A stands for all integer linear combinations of 
the fundamental weights of su(N). The remaining integrations are immediate, the 
integral over dA^ 1 producing S((f/ l + A), so that indeed only constant (jJ 1 configurations 
contribute. One can now use the classic formulas [12|, [TJ], |I5| , 

( ic i tun ' '' 1 1 " ' ' ^ 

and 



detM^W-nCM h ) (15) 

<5 



C 2 (/i) = (// + p) • (/i + p) -P-P (17) 

where dim(/i) and C*2(p) is the dimension and quadratic Casimir of the irreducible 
representation of SU(N) with highest weight //, 5 (<5 + ) is the set of all (positive) roots, 
and p is the half-sum of the positive roots.0 



Ji 



As remarked in [11, 12] the Weyl determinant in (15) vanishes whenever (j> lies 



in the walls of a Weyl chamber. Thus the integration should be restricted to regular 



3 The quantity (/x + p) coincides with the variable a defined by Okubo p/7j j , his eqn.(15), in his 
study of third order Casimir operators. 



valued 0^ fields, i.e. those in the interior of the Weyl chambers (see Jl3| for a detailed 
discussion). In the source free theory all (fJ 1 can then be rotated to the interior of 
the fundamental Weyl chamber, and the factor gets cancelled. Every element of 
A then corresponds to a unique element (p + p) in the interior of the fundamental 
Weyl chamber. However when the source J is not zero, one must retain the average 
over the Weyl group. In fact, that is why we have kept the harmless constant of 
the Weyl integral formula in front of the path integral (T2). The result for the su(N) 
sector can then be written, 

A* I I ere IV 

(18) 

where the sum goes over representations of su(N) compatible with u(l) charge n as 
explained above. In ( |i8| ) we have fixed the allowed renormalization terms 0, [12|] of the 
form exp(coe 2 A) and exp(ci(2 — 2g)), to give agreement with the usual normalization 
chosen for J = 0. We can now combine the u(l) and su(N) pieces to give the total 
partition function for U(N), 

Zx g (J) =Edim(/) 2 - 2 ^exp(^^l) {exp^ dp Tr( J) 2 ) 

>W E exp(*e 2 / E9 ^(a(Z + p), J*))} (19) 

where now I runs over the highest weights of irreducible representations of U(N). 
(The V are the row lengths of the U(N) Young tableau for the representation I.) It 
should be emphasized that the last term in (|T9"D is a consequence of the coupling of 
the source to the nontrivial T topological sectors of the gauge field. Notice that for 
U(N), the Weyl group W is the symmetric group Sjy. This computation can also be 



trivially extended to the case of the generalized 2d QCD of |18|, ^5|. In fact, since 
we have seen that only constant configurations of t^ 1 contribute, our calculation also 
applies when higher powers of 4> are present. 
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3 The Correlators 



In this section we will study the field strength correlation functions, obtained from 
functional derivatives with respect to J in (J). Since the average over the Weyl 
group plays an essential role, it is convenient to use a basis for h@u(l) where the 
action of the Weyl group is as simple as possible. For k we will use the usual generators 
of the fundamental representation of SU(N), and for the generators of T we will use 
the diagonal matrices En i = 1,...,N, where En has only one nonzero entry, 1 on 
the i th position of the diagonal^. We adopt the normalization Tr (T a T b ) = 5 ab . The 
elements of the Weyl group Sn then just permute the entries of the diagonal matrices 
in h © u(l). (The Weyl group acts as the identity on u(l).) 

As expected from the symmetry £ — > — £ of ([I]) with J = 0, all odd point 
functions vanish. In our framework, the odd point functions vanish since various 
products with an odd number of factors of the form, 

(i+pr d+pv 

are averaged over W . However, for each represention with highest weight I there 
exists a conjugate representation with highest weight I such that (l+p) = —a(l + p), 
where Oi — N — i + 1. Moreover the dimensions and quadratic casimirs of I and I are 
the same, so that the contribution of / cancels that of I in the sum over I in the odd 
point functions, and therefore these vanish. (For self-conjugate representations one 
will have (l + p) = —a(l + p) so that the average over the Weyl group also gives zero). 

There are two terms which depend on J in Zs g (J). The term exp[/ Sg dpTr(J 2 )] 

gives rise to contact terms in the expectation values of the field strengths. If one 

wants to consider products of fields at the same point one needs to regulate these 

terms. This term is closely related to the arbitrary renormalization term exp(coe 2 A) 
4 In this basis one has that the l l are the row lengths of the Young tableau of the U(N) represen- 
tation I and that p l = — i. 
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that appears in front of the partition function. In fact by differentiating the J = 
partition function with respect to e 2 , we can conclude that the normalized correlator 

{^ g d^e{x)) = 2e 2 A^ x F (20) 

where F = log Zj^ is the free energy.0 Therefore in order to bring fields to 
the same point we must normalize the contact term accordingly. This is done by 
adjusting the coefficient Co in the renormalization term. 

For example, with this renormalization prescription the regularized 2-point func- 
tion becomes 

(Cix^y)} = ^Edim(/) 2 -%xp(^^) 



Z s 9 , 



X 



(p,p)S ab x2 1 



x , ^ y - Wl Ui+pr a (i+pr b ] (2i) 

where 5l y is 1 if x — y, and zero otherwise. Notice the (I + p) terms are present only 
when a and b are h © tt(l) indices, i.e. when the indices lie in the Cartan subalgebra. 
The average over the Weyl group can be done explicitly with the result (see appendix 
B) 

1 + py a ( i + py b = p ab ( i + py + mabn2 ( 22 ) 



\w 



where 



and 



p a b={ N(N-l) [ia ^ b ^ m a b= \ N(i-1) iia ^ b 

if a = b j if a = & 



n = y~] V = total number of boxes in the tableau defined by I. 

i 

Thus we have the final result for the regularized 2-point function 



5 All correlators will be normalized by dividing by the partition function. 
That is {fix) ■ ■ ■ e(y)) = z^o)oj^--- eA^)Z Sg (J)\ J=0 
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(ewem = ^y.^) 2 ~ 29 ^m eM 2 Ca(0 ) 

^9 I " 

xl^0^%,y-(P ab (l + PT + ™ ab n 2 )} (23) 

This is finite as x — > y. From here ( p0|) follows as well. Notice that the correlator is 
essentially topological in the gauge we have chosen. 

Similarly the regularized 4-point function is, 

(?(x)?(v)?(*)t») = ^Edim(/) 2 - 29 exp(^^l) x 

^9 I 

{[MjW. + 2 permut.)] - E + P)°*(! + pT + 5 P ermut "] 

(once again the (I + p) terms are present only for h(Bu(l) valued indices in the Cartan 
subalgebra). The average over W gives (see appendix B), 

1 



Y^(i+ P r(i+pp(i+pr(i+ P y 



a ijki ^(l + p )4 + b W[(i + p )2]2 + c m n J2(l + pf + ^' fcZ n 2 (Z + pf + e^'n 4 (25) 

i i 

where the coefficents a lJ , etc can be found in equation (B.6). Combining (p3| ) and 
f ^5|) gives the regularized 4-point function 

,8 



?(x)?(v)t h (z)?(v>)) = ^ 1 Edim(/) 2 - 2g exp( eMC M - 



+ 2 permut.)] - ^^[^ >y (p«(Z + p) 2 + m A V) + 5 permut.] 
+ [a l]kl Y J {l + P)t + b ijkl [(l + pf] 2 + c l3kl nJ2(l + P?i + d ljkl n 2 (l + pf + e^ M n 4 ]}(26) 

i i 

Once again we see that the correlator is topological. Contracting with the group 
generators produces a gauge invariant quantity. (Some useful identities for traces and 
higher order Casimirs can be found in [HJ). We have 
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<Tr£ 4 (*)> = ^Edim(/) 2 - 2 *exp( 1^1 ) x 

^9 I 

— ( M 2 — 1 ") 2 1 /V 2 — 1 /V 2 — 1 1 

[_^! iL(Ar 2 + 6N + -) - -(TV - l)C 2 m - ^n 2 + ±C 4 (Z)1 (27) 

where C±(l) is a fourth order Casimir for U(N) defined by C±(l) = YUW + P)*] 4 - 

The structure of the higher point functions is now clear. For the 2p point func- 
tion there's a Casimir invariant C*2 P , and products of (p, p) with lower order casimirs 
of even order. The gauge invariant quantities (Tr£ 2p (x)) have no dependence on x 
at all, while the non-invariant correlators are also topological and only depend on 
contact terms in a trivial way. Thus, this gauge could be called a topological gauge. 
In fact the theory is invariant under diffeomorphisms of S 9 that preserve the measure 
djji,. Thus the x dependence of gauge invariant quantities could come only from factors 



of J g(x) which were already absorbed in the definition of £ in ([!]). Similarly, in the 
computation of Wilson loops only the area enclosed by the loops was relevant (see 
for example [l(|). However, notice that even the non-gauge invariant uncontracted 
correlators (plD,(p4|) above have only a trivial position dependence coming from the 
contact term. Therefore the gauge (0) is particularly well suited to exhibit the topo- 
logical nature of the theory, in the sense that expectation values of local non-gauge 
invariant quantities in this gauge are also independent of position. 



4 The Phase Transition On S* 2 

The partition function on the sphere is, 

Z = £dim(i?) 2 exp(^^) (28) 

R 

where now A = X(Area), with A = e 2 N held fixed when iV — > oo. Although the 
exponential decreases rapidly for large representations, for A sufficiently small the 
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series diverges due to the dimension term raised to a positive power. In fact, Douglas 
and Kazakov [I5| have shown that in the large N limit the theory undergoes a third- 



order phase transition at A cr = tc 2 . This phase transition was also explained in 
pOfl as a consequence of the presence of instantons (classical solutions of the field 



equations). When A approaches it 2 , the contribution of the instanton configurations 
to the partition function becomes dominant, and the phase transition is induced. 
The relation to the string interpretation was also explored in |2^, |23|. It has been 
shown that a similar phase transition takes place for YM2 on the cylinder pi |, on the 
projective plane RP 2 , and for Sp(N) and SO(N) on the sphere as well, p4|. Further, 
the generalized 2d Yang-Mills theory possesses a rich phase structure, with phase 
transitions possible for g > 1, and for particular regions of the coupling constants 



space p5f . The Wilson loops in the large N limit have been studied in both the weak 



and strong coupling phases |L9L E71 ER 29, 24 



Here we will study the phase transition on the sphere for the electric field 
correlators of YM 2 in the large N limit. Let us briefly review the methods and known 
results. The irreducible representations of U(N) are labelled by the Young tableau 
row lengths I 1 >l 2 > .... > l N . Also, 

(l + p y = _±i -i + V , i = 1...N (29) 



In the large N limit, these variables can be replaced by continum variables [26 



1 

x = — with < x < 1 
N ~ ~ 

K x ) = jj and Y. = N Jl dx ( 30 ) 



It is then useful to define, 



h(x) = -l(x) + x - i (31) 



When one expresses the partition function (28) in terms of these continum variables, 
the sum over representations becomes a quantum mechanical functional integral for, 

Seff(h) = - Cdx I" dy \og[h(x) - h(y)\ + £ /* dxh 2 (x) - ^ (32) 
Jo Jo 2 Jo 24 
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with 



Jvhexp(-N 2 S eff (h)) (33) 



Since N is large, this can be computed by the saddle-point approximation |26], |T9 
Define the density of boxes in a Young tableau by 



= where u (h\ < i (34) 

oh 



The saddle point equation is 



P / d s^- = ^ (35) 



For A < A cr = 7r , the solution to this equation is given by the semicircle law p6| , 



U we ak{h) = ^J- h2 > With < k < VI ( 36 ) 

However for A > A cr = tc 2 the above solution violates the condition u(h) < 1 in 
and a different solution must be sought [19]. The ansatz 



gives the solution 



, 1 if -b < h < b 

U s trong(h) = ^ (37) 

u(h) otherwise 



f(h) = I ds^ 



-h- - log^-4 + J (a 2 - h 2 )(b 2 - h 2 ) I' ; dS = (31 

2 h + b V J-b (/ l _ s )^( Q 2_ s 2)( & 2 _ s 2) 



Following |l9l we will now present the behaviour of the free energy, F = -^-logZ, at 
the critical point. This will enable us to review the method and notation needed to 
describe the phase transition of the correlators. We have, 



dA OA 2 Jo 24 
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where one uses the fact that u(h) is an even function both in the strong and weak 
coupling solutions. When A < ir 2 from 



F\A) 



1 

24 



1 

2A 



(40) 



For the strong coupling phase we need to use the large h expansion of (|38[) . The 
0(h) and 0(h~ r ) terms give, 

a 



AK 
~A~ 



and 



(41) 
(42) 



A = 8EK - 4k 2 K 

where k = b/a is the modulus of the elliptic integrals, k! = \/l — k 2 is the comple- 
mentary modulus and K = K(b/a), E = E(b/a) are the complete elliptic integrals 
of the first and second kind respectively (see for example ||30|| ). To explore the be- 
haviour near the critical point, it is useful to express these quantities in terms of 



theta functions on a torus of complex modulus r [IE, 3C[. One sets q = expiirr and 
6i = 9i(0\r), i = 0, 1,2,3. When we approach the critical point, b — > and q — > 0. 
We then define 5 = A — A cr = A — it 2 . Expanding (|42|) in q gives, 

A = 7r 2 (l + 8q-8q 2 + 32q 3 

and inverting for 5 we find, 

5 ^ 5 



(43) 



Q 



8tt 2 



1 + 



+ •••) 



From the weak and strong coupling expansions of (f40|) one gets 

1 (A-A c ^ 2 



strong 



weak 



5 2 

7T 6 



+ ■ ■ 



7T- 



7T' J 



(44) 



(45) 



which is the third-order phase-transition for the free energy of [1^] . 



From (|20D , section 3, in the large N limitQ (note that here F = -™ \ogZ) 
-(Tre 2 (^))=2A— = 2AF'(A) 



(46) 



3 In this section 2 will denote a point on S 2 and x will always be the continuous index x 



N ■ 
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and therefore the gauge- invariant two-point function behaves as F' (A) at the critical 
point, with a second-order phase-transition. 



The gauge invariant four-point function was given in fl27|). In the large N 
limit, only the critical representation survives at the saddle point. Moreover, since 
the critical representation is self-conjugate, h(x) is odd about x — ~, from fl3T|) the 
corresponding u(l) charge n = f l(x) vanishes. Notice that, 

{I + P)(x) = + P)' = -Kx) 



then 

N 



(T< 4 W)=A 4 {— --[J Q dxh 2 (x)--] + l j o dxh\x)} (47) 
In the weak phase we can use ( |36D to obtain, 



1 2 

dxh weak (x) = — (Ah 



From the 0(h ) term in (|38|) in the strong coupling phase we have, 



[ dxh% rong (x) = f^(ef - eje! - « + ef) + g£(6>§ + \e\e\ + el) (49) 



Expanding for small q and using (f44|) we find 



dxh str(mg {x) - / dxh weak (x) = — ^— + ■■■ (50) 



so that the phase-transition for the four-point function is also second-order. 
Similarly we can use the 0(h~ 7 ) term in fl38|) to show that, 

-/ 

Jo 

W [ 2688^ 3 + ° 2 ] ~ 672^ 3 ° 2 + ° 2 ° a) " 448^ 2 " 48 Wa " 96 W + ^ )] 

(27r) 6 4 12 12 5 m8 „4 /.SMx , 7 z)8/z)4 , z)4> 



d> x h strong (x) 



-[-^3 + * 2 ") + ttt(« + m) + izFM + m (51) 



^ 6 l 2i ^ 1 ' 42 v d 2 2 dy 48 
so that, 

dxh% (x) - f dxh 6 weak (x) = — ^— + ■ • • (52) 

JO 7T° 7T 4 
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Therefore the six-point function also has a second-order phase transition at the critical 
point. (In fact it has been shown that 



1 ^ rl ^ 

d%h strong (x) — dxh weak {x) 



has a second-order transition for any n, so that all higher-point correlators of the field 



strength have second-order phase-transitions as well |31|) 



The master field for the field strength and gauge potential, on the sphere in the 
large N limit, based on the above results, which reproduces the regulated correlators 



is presented in [32|. The unregulated correlators, that is including non- regularized 



contact terms, can also be described by an explicit master field representation [36 



5 Application to Wilson Loops 

One of the motivations for studying field strength correlators is to understand the 
couplings of matter to 2D Yang-Mills theory on Riemann surfaces. To illustrate this 
direction of research, we compute expectation values for contractible Wilson loops 
using a version of Stokes' theorem, equation (|53|) , and our correlators and compare 
these calculations with known exact results. Although we do not have a complete 
proof of the validity of this version of the theorem, our results are so compelling, 
that we include them so as to stimulate further research. To be specific, we compute 
contractible Wilson loops on T, g for U(l) and SU(2) gauge theories, using Stokes' 
theorem to relate the holonomy of the gauge field in the representation r around a 
closed loop to the exponential of the integral of the gauge curvature over the interior 
of the loop. That is we consider 

(Wrfo)) = (Tr Vexp(i I A)) = (Trexp(i / F)) (53) 

Jf Jd 

where we will take a contractible, non-self-intersecting loop 7 = 3D. Obviously ( |53"D 
is valid for a U(l) gauge theory. However, motivated by the abelianization of the 
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partition function, we attempt to use (|53|) for the SU (2) gauge theory as well. As we 
describe below, we have computed (^(7)) for SU (2) up to 0(A 4 ), where A is the area 
enclosed by the loop 7, and verified that the Stokes' theorem (|5^) is valid for this case 
as well, in the Blau-Thompson abelianization gauge. It should be emphasized that 
(p)5|) appears to work for non-abelian theories, with gauge group SU(2), computed 
in the abelianization gauge, as a result of explicit calculations. Unfortunately, we 
cannot provide a proof of this remarkable result, so that the application of (^) to 
SU(2) non-abelian theory in this gauge remains a conjecture. We note that for non- 
abelian Yang-Mills theories in general gauges, one must use a non-abelian version of 
the Stokes' theorem, where one of the variables in the surface integral is ordered and 
the gauge curvature appears conjugated by holonomies, so that gauge invariance is 
preserved p3 |. By contrast, in (|53]) there is no ordering in the integration over the 



disc and F is not conjugated by holonomies. Nevertheless, it apparently computes 
the Wilson loop in the abelianization gauge. 

The expectation values of Wilson loops in two dimensional Yang-Mills theory 
are known; see for example |l(J. For a loop 7 = 3D on E fl , in the representation r 



we have the exact result 

W , ~e 2 AC 2 (Z) ^ e 2 A(C 2 (l) - C 2 (A)) 
{W r {j)) = 2^dsm{l) 9 exp( ) ^ dim(A) exp[ ] 

l A \el®r A 

(54) 

where A is the area of the disc D enclosed by the Wilson loop. The exponential in 
([54]) is expanded in powers of A, which then expresses (W r {l)) as a power series in 
A. On the other hand, one can expand the right-hand side of (^) as 

(W r ( 7 )> = Tr(l) - i(Tr ^ F(x) ^ F(y)} + i(Tr( ^ F) 4 ) + • • • (55) 



where we have used our result that odd correlators of field strengths vanish. 
Comparing the expansion of (|54j) in powers of A with ( |55"D seems to present a para- 
dox, as (0) contains odd powers of A, while (|55|) only involves an even number of 
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integrations on the disc D. In fact, there is a very elegant resolution of this issue 
when one uses the (unregulated) correlation functions of section 3, as these corre- 
lators have contributions from the contact terms, obtained by means of functional 
derivatives with respect to the source from the Tr(J 2 ) term in ([191) . We will refer to 
contributions from this term, as the contact terms of the correlators, and from the 



linear term in J in (19) as the topological terms. 



For example, the contact term (~ J g(x) J g(y)S 2 (x — y)) in the two point func- 
tion will produce a linear term in A in fl55"D when integrated over the disc. This is also 
the case for the Wilson loop on the plane. See for example Bralic in In general 



each contact term and each topological term in a field strength correlator produces a 
term linear in A when inserted in (1551). 



We now turn to the explicit verification of fl55|), correct to 0(A 4 ), with U(l) 
and SU(2) considered separately for clarity. 



5.1 U(l) gauge theory 

The Wilson loop in an arbitrary representation r of U(l) can be written as 

nxrt w l -e 2 An 2 -e 2 A(r 2 + 2nr) 

(W r (j)) = £ exp( ) exp( K — '-) (56) 

where n runs over all integers. 

The expansion of (|56"D in powers of A yields 



00 p 2 An 2 

(W r ( 7 )> = E exp( — ) • 

n=— oo 

.11 _ f! r 2 A + If!( r 4 + 4n 2 r2) A 2 _ 1 ^,6 + 12n V)A 3 + 

1 2 2! 4 v 1 3! 8 V ; 



+^(r 8 + 24n 2 r 6 + 16n 4 r 4 )A 4 + ■ • ■} (57) 
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= (1) - -A(r 2 > + -A 2 (r 4 + 4nV) - — A 3 (r 6 + 12nV) + 
2 8 48 

+ I^A 4 (r 8 + 24n 2 r 6 + 16nV) + --- (58) 
4! 16 X l \ ) 

where the expectation values in ( |58"| ) are defined in the obvious way from fl5"T|). 

Our objective is to compute the right-hand side of d55|) and to compare it with 
(j58|). The relevant generating functional of correlation functions is (|6]), specialized to 
N — 1, and with each F multiplied by r in accord with our normalization conventions. 
Let us sketch the explicit comparisons, using (]55| ) a nd our unregulated correlators. 
For convenience we write F c or F t for the part of the field strength correlators coming 
from the contact terms or topological terms in the generating functional. We indicate 
these terms in the expectation values schematically, as these usually correspond to a 
number of permutations of contact and topological terms. We stress that F c and F t 
have only symbolic meaning, as labels of terms in the correlation functions. 

0(A) : 

A f dx f dy (F c (x)F c (y)) = (r 2 )(-i) / d^x) f d/x(y) e 2 5 2 (x - y) = -^(r 2 )A 

2 JD JD 2 JD JD 2 

(59) 

which agrees with ([58]). 
0(A 2 ) : 

-\l dx I dy (F t (x)F t (y)) = (*e 2 ) 2 r 2 (-W>A 2 (60) 

2 JD JD 2 

and 

l { / D ) 4 ( F c) = i[ < r V( J/[S 2 (x - y)5\z -w) + 2 permut.] = ^e 4 (r 4 )A 2 (61) 

The sum of ( |60"D and ( |6T| ) agrees with the 0(A 2 ) term in (|58|). Further, the two terms 
in (|58|) of order 0(A 2 ) originate from fl6lD and (|60| ) respectively. 



0(A 3 ): 



H ( / D )4(Fc2F * 2) = I[ r4 (- e4 )^ 2 )( / D ) 4 ^ 2 [5 2 (a:-y) + 5 permut.] = -^r 4 (n 2 )A 3 (62) 
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and 



-1( / ) 6 (F C 6 ) = -^(r 6 )( ( f e 6 [5 2 (x-y)6 2 (z-w)5 2 (v-u)+U permut. 
6! Jd 6! Jd 



Again, the agreement with fl58|) is evident. 



- — (r 6 )A 3 
48 V 1 

(63) 



0(A 4 ): 



while 



and 



47 ( //^ 4) = 4T r4 ^ )A4 



-A I ) 6 (F 4 F t 2 ) = -r 6 (n 2 )A 4 

6rV w * 7 16 v 7 



4!16 



[r 8 )A 4 



Thus, (|64T) to (166|) agree with the 0(A 4 ) term of (|58|). 



(64) 



(65) 



(66) 



Therefore, we have shown how to obtain the expansion of the Wilson loop 
expectation value in powers of the area, for a contractible non-self-intersecting loop, 
directly from the field strength correlation functions. 

We note that each term in (|58|) comes from a distinct combination of contact 
and topological terms in the expectation values. This sets the stage for a similar 
calculation for SU(2). 



5.2 SU(2): Wilson Loop in the Fundamental Representation 

The Wilson loop expectation value of a gauge field in the fundamental representation 
for a contractible non self-intersecting loop is 

w(t)> = 2^ dim U) ex P( 9 ) 2^ dim(A)exp( ) 

(67) 
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where dim(j) =j + l, C 2 {j) = + 2) and j®f = {j — l)@ (j + 1) for j > 1. The 
expansion of ( p7| ) in powers of A gives, after some elementary but tedious algebra 

(W>( 7 )> = ^dim(j) 2 - 29 exp(-^Cz)) • 

■{Tr(l) - 3 -e 2 A + ^[4(j + l) 2 + 5]A 2 - ~[2D{j + I) 2 + 7]A 3 + 

+ ^^[16(J + I) 4 + 56(j + l) 2 + 9]A 4 + ■ ■ -}(68) 
= (2) - ^e 2 A(l) + ^A 2 (4(j + l) 2 + 5) - ~^A 3 (20(j + l) 2 + 7) + 

+ ^A 4 (16(j + I) 4 + 56(j + l) 2 + 9) + ■ • • (69) 



The generating functional for the correlators of field strengths in the fundamen- 



tal representation is given by (18). For 377(2) it is convenient to carry out the Weyl 



sum in fll8"|) directly, with the result 

= E dim(j) 2 - 2 ^ exp( 2 ° 2Uj ) exp(^-/ Sg ^Tr( J 2 )) cos(^=(j+l)/ Ss dptJ ) 

(70) 

where ^^'c 3 is the diagonal component of J, with our choice of normalization. 
Let us compute the right-hand side of ( [55)) using the correlators generated by (fT0|). 

0(A) : 

-\fdxf dy(Tr(F c (x)F c (y))) = e 2 (5 ab Tr(T a T b ))(-l-) f d»{x) f d^y)5\x-y) = -^e 2 (l) 
2 Jd Jd 2 Jd Jd 2 

(71) 

which agrees with QB5|). 
0(A 2 ) : 

(^) 4 (Tr(F c 4 )) is obtained from e A [5 ij S kl 5 2 (x - y)5 2 (z - w) + 2 permut.] 
which must be contracted with traces of the generators with structure Tr (T l T J T l T J ) 
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and (Tr(T*T*)) 2 to obtain 

^( / D ) 4 (Tr(F c 4 )) = ^e 4 (l)A 2 (72) 

Similarly 

- ~( //<Tr(F 2 )> = ^((j + l) 2 )Tr[-La 3 i=a 3 ]A 2 = ^((j + 1) 2 )A 2 (73) 

Equations ([72]) and (|73|) agree with the 0(A 2 ) term in (|69D. (Note that only generators 
in the Lie algebra of the maximal torus contribute to the topological terms.) 

Q(A 3 ) : 

The calculations become increasingly lengthy, so we only summarize the results 
I\ ( J/^ F " F ^ = " e6 ^2 A3((J + 1)2) (?4) 

and 

-^(//wo = -Ilr 8 a>A 3 < 75 > 

Equations ( |74"D and (|75|) agree with the 0(A 3 ) term of (|69|). 
Q(A 4 ) : 

^(//(Tr(F/)) = ^A 4 ((j + l) 4 ) (76) 

while 

and 

-, ( / ) 8 (Tr(F 8 )) = -i— e 8 (l)A 4 (78) 
8!W X V c " 41128 W V ; 

The computation of ([76]) to ( [78] ) is very lengthy, but nevertheless agrees with 
(|69|). Note that in the expectation values F 2 gives e 4 (j + 1) 2 A 2 and F 2 gives e 2 A up 
to an overall constant in the computations ([H]) to (|7|). This clearly generalizes to 
arbitrary correlators. 
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In summary, we have verified that the contractable non-self-intesecting Wilson 
loop for U(l) and SU(2) can be computed from our correlators and the right-hand 
side of (|55|), at least to 0(A 4 ). One would think that a discrepancy would have 
already appeared at or before 0(A 4 ), so that this provides an impressive check of our 
correlators, and also "experimental" evidence that the Stokes' theorem ( pB| ) is valid 
for the abelianization gauge, even for SU (2)Q A formal proof of this remarkable 
result is lacking. Finally, we mention that we have not assumed that the loop 7 
was infinitesimal, or that A was small. This provides one more advantage of the 
abelianization gauge, which seems to provide an enormous simplification as compared 
to the non-abelian Stokes' theorem f33fl . 
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Appendix A: The Kernel on the Disc 

In the main text we started with the action for 2d Y.M. coupled to an external 



source on the closed surface S 9 , and then used the abelianization method of JTT], [12 



for the path integral, where the nontrivial U(N) and T bundles over H g had to be 
incorporated. Functional differentiation with respect to J was then used to derive the 



correlators. Similar results can be obtained by a different method. In |T0[ Blau and 
Thompson computed the kernel for 2d Y.M. on a disc and used it to find the partition 
function and Wilson loops for closed surfaces. One starts with the theory on a disc 
and sets a boundary condition for the holonomy of the gauge field around the bound- 
ary 7(t) of the disc D. The path-ordered exponential representing this holonomy is 
expressed in terms of an ordinary exponential, by means of a functional integration 



7 For SU(N) with N > 2 the abelian Stokes theorem works up to order A 2 but seems to fail at 
next order E}3]. 
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over auxiliary anticommuting fields. This introduces an explicit dependence of 
the integrand in the path integral. However, for the simple topology of the disc, the 
Schwinger-Fock gauge is available, and this allows one to express in terms of F^, 
so that the Nicolai map becomes useful. The integration over the gauge field and over 
the auxiliary fermionic fields can then be performed |I0 . 



When we couple the theory to an external source J, the calculation follows 
along similar lines. The resulting partition function is, 

Z D (g, J) = £ d im(r){exp(^^) exp ^ dp Tr( J 2 )] XrfrPt exp^ f p£r a ))} 

7 (A.l) 

where r denotes an irreducible representation of the gauge group G, Xi an d r a are 
respectively the corresponding Weyl character and generators, g is the holonomy 
around j(t), Pt is the path-ordering operator and p° is given by 

p£(t) = f sdse^(t)JgJ^)J\ S1 (t)) (A.2) 



where we have cartesian coordinates x 1 ^ 2 on the disc D with, 

a?(s,t) = srf(t) for 0<s,t<l (A.3) 

When glueing together two such discs, by identifying holonomies along bound- 
aries and integrating over all such possible holonomies, one should be careful in 
specifying the t — and t — 1 points of the boundaries consistently. Results for 
closed surfaces are obtained by decomposing the holonomy into the product of the 
holonomies along the cycles of S 9 and making the corresponding identifications. One 



can now perform functional derivatives of (|A. 1|) with respect to J to compute the field 



strength correlation functions. The correlators should be independent of the way one 
chooses to decompose the surface S s . So, for example, one can construct a sphere 
out of two discs, with boundaries j(t) and 7 _1 (t), or by considering only one disc 
with boundary aa~ l {t), and the answer should be the same in the two cases. Let 
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us construct the regularized 2-point function on the sphere, considering the 1 disc 
approach first. One decomposes the holonomy in to gg~ 1 in ( |A.l ) and in this case 



the integration over holonomies with J G dg just produces a factor of unity. Functional 
differentiation with respect to J then produces, 

(Cix^iy)} = PSlj ab + ^(- e4 ) E dim W ex P( ^ Ai )Xr(r a r 6 ) (A.4) 

where (3 is a renormalization dependent constant. Notice this is a different gauge 
than the one considered in the body of the text. If we construct the sphere out of 
two discs the partition function will be, 

Z s z= dg Z D (g,J D ) Z jy {g~ x , Jjy) (A.5) 

J G 

and this integration can be carried out using the familiar orthogonality relations 
between Weyl characters, 

dg Xr(gh)x T >(g- l h') = -^-Xr(Wi') (A.6) 
G dim(rj 

with the result, 

z s 2 =E dim ( r ) ex P( ^ — ) exp (y J s2 dfiTr(J 2 )) x r (Ptexp(ie 2 ^ p)P t exp(ie 2 j ; p)) 

(A.?)" 

where p M = Yl.aP^ ' an d p is obtained from J D > as in ( |A.2| ). The two-point function 
obtained from ( |A.7| ) consists of several terms which depend on wheter or not x and y 
are in D or D . These terms can be understood as arising from glueing the one-point 
function on D with the one-point function on D' , and the two-point function on D 
with the partition function on D and vice-versa. The result is the same as the one 
in ( |A.4| ), as expected. Notice that in this calculation, the t valued fields play no 
special role, all indices in the Lie algebra being treated equally. The gauge invariant 
two-point function, for the appropriate choice of the renormaliza tion constant /?, is 
identical to the one computed in section 3. 
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For higher genus surfaces and higher point functions this method becomes cum- 
bersome. The explicit dependence of the partition function on the choice of coor- 
dinates s and t appears in the correlators, the ordering in which the generators r ai 
appear inside traces being dependent on the ordering of the coordinates t(xi). Thus 
the symmetry must be restored by hand, and only then do the higher order Casimir 
operators appear. So, for closed surfaces the method based on abelianization is clearly 
simpler and more elegant, however this second method can be applied to the compu- 
tation of electric field correlators on surfaces with boundaries. 

Appendix B: The Weyl Group Averages 

Here we will present the averages over the Weyl group used in section 3. Recall 
that the Weyl group acts as the symmetric group of N elements, Sjv, by permuting 
the coordinates (7 + p)\ All the identities that are needed can be deduced from, 

JV 

£(/ + p)' = n (B.l) 

i=i 

where n is the u(l) highest weight corresponding to the U(N) irreducible represen- 
tation of highest weight I. For the two-point function ([Hp we need to compute, 



If a 7^ b (p.2| ) becomes, 



and then the sum over Oh produces 

1 



so that finally we have for the case a ^ b 
1 



_5>(z + p)*.-((z + p)*.; 



£(i + P y°(i + P y» = AUAT A n 2 - (i + pf] 



\W\ V N(N-l) 



28 



(Note we have the normalization v 2 = J2i=o vlv \ f° r v £ £)• 



Similar reasoning for a = b gives the result presented in section 3: 
1 



with, 



P 



ab 



w 



+ PY a {l + PT" = P a \l + P? + ™ ab n 2 



n(n-i) if a 
it if a = & 



and 



m 



, ab 



n(n-i) if a b 
if a = b 



(B-3) 



(B.4) 



For the four-point function we have by similar arguments, 
1 



W 



ijkl 



E( Z + P)i + bm [( l + P?f + ci3Un E( Z + P)l + d ijkl n 2 (l + p) 2 + e^'n 4 (B.5) 
v p 

The coefhcents of the various Weyl group invariant terms are completely symmetric 
in the indices i,j,k,l. With 

1 i = j 
i^j 

we have (with no sum on repeated indices in flB.6|) ) : 



ijkl 



-6(AT-4)! 
N\ 

2(N-3)\ 

N\ 
-(N-2)\ 

N\ 
-(N-2)\ 

N\ 

m 



i,j,k,l all different 

S ij £ ik £ u £ ki 

5 ij 5 jk e il 
5 ij 5 ki £ ik 

i = j = k = I 



b ij 



ki 



Jjkl 



8(AT-4)! 
AT! 



i,j,k,l all different 

A -■"■ii fiij £ ik £ il £ kl 

otherwise 



JV! 
(AT-2)! 



A7 



JV! 





3(A/-4)! 
N\ 

-(N-3)\ 
N\ 



(N-2)\ 
AT! 



-6(jV-4)! 
AT! 

(jV-3)! 
AT! 





i,j,k,l all different 

S ij e ik £ il e kl 

5 ij S jk e il 
5 ij 5 ki £ ik 

i = j = k = I 

i,j,k,l all different 

5 ij £ ik £ u £ ki 

otherwise 



AT! 



i,j,k,l all different 
otherwise 



(B.6) 
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Notice that when one contracts with the generators of the maximal torus in (|27|), 
only the i = j = k = I case contributes. 
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